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Abstract. Annihilation radiation from supersymmetric particles at the Galactic 
center could be greatly enhanced if the dark matter density is peaked around the 
supermassive black hole. Arguments for and against the existence of density spikes are 
reviewed. Spikes are destroyed during mergers, and there is strong evidence for this 
effect in stellar density profiles. The dark matter spike at the Galactic center probably 
suffered this fate. 

Indirect detection schemes are based on searches for gamma rays, neutrinos or other 
annihilation by-products from supersymmetric particles in the dark matter halo of the 
Milky Way [1]. The flux depends on the squared density of particles integrated along 
the line of sight and the signal is greatly enhanced in directions where the dark matter 
is clumped. This includes the center of the Milky Way where the density in a smooth 
halo would be maximum[2, 3]. The signal from the Galactic center is further enhanced 
if there is a dark matter "spike" associated with the central supermassive black hole 
(SBH)[4]. Growth of a SBH causes the orbits of nearby dark matter particles to shrink 
and their density to rise; for a wide range of initial conditions, the resulting profile is a 
steep power law, p ~ r~ 2 , implying a formally divergent annihilation flux from near the 
SBH [4]. 

This article reviews the arguments for and against the existence of dark matter 
density spikes. A useful guide to the distribution of dark matter at the very centers of 
galaxies is the stellar distribution; in the absence of gaseous dissipation, stars act like a 
nearly collisionless fluid and should react to the presence of a SBH in the same way as the 
dark matter. Steep spikes in the stellar luminosity profile are only seen in galaxies which 
show evidence of dissipative formation. Other galaxies exhibit shallow inner profiles, 
and this is consistent with a model in which binary SBHs injected energy into the stellar 
fluid during the mergers that formed the galaxies. Mergers almost certainly occurred 
during the formation of SBHs and their host bulges, implying destruction of the dark 
matter spikes. 



1. The Adiabatic Growth Model 



If a black hole grows at the center of a collisionless fluid, the density around it also 
grows as the black hole's gravity causes the surrounding orbits to shrink. The result is 
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Figure 1. Spike formation by adiabatic growth of a black hole, (a) Thin curve is the 
initial model, formed by collapse in an iV-body simulation. Heavy curves are density 
profiles after growth of a central "black hole" containing 0.3%, 1% and 3% of the total 
mass. The reference line has a logarithmic slope of —2. (From Ref. 5.) (b) Initial 
models were spherical isotropic halos with power-law density profiles, p oc r~ 70 ; 70 
increases upwards in steps of 0.25. The radial scale is normalized to in the initial 
halo. The slope of the final profile at r < r/j is almost independent of the initial slope. 



simplest to compute under the "adiabatic" approximation in which the growth time is 
long compared with orbital periods; this is reasonable in the case of SBHs at galactic 
centers, whose growth time is thought to be ~ 10 8 yr compared with orbital periods 
of ~ 10 6 yr. The resulting density profile depends somewhat on the initial state. One 
possibility is a constant-density, isothermal core; the resulting spike is a power law, 
p oc r~ L5 , r < rh = GM,/a 2 , with M. the final black hole mass and a the ID particle 
velocity dispersion before the appearance of the black hole [5, 6]. 

A more reasonable guess for the initial profile is a power law, p oc r~ 7() ; power-law 
central profiles are generic outcomes of hierarchical structure formation simulations [8, 9], 
although the index of the power law is debated. Adiabatic growth in a power law gives 
a final profile that is also a power law near the black hole but with steeper index[4]: 

P/(r)cxr-T, 7 = 2 + -^- (0 < 70 < 2). (1) 

4-70 

The final index is almost independent of the initial index: 2.25 < 7 < 2.5 (Figure 1). 
Other reasonable initial profiles can be found[10] that generate final profiles which fill 
the gap between the power laws and the isothermal sphere. Hence in the absence of 
detailed knowledge about the initial state, the adiabatic growth model predicts density 
profiles such that 

p f {r)^r-\ 1.5< 7 <2.5. (2) 
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Even slopes at the low end of this range should not be ruled out since the latest 
simulations [11] suggest that 70 may decline monotonically toward the center. Essentially 
nothing is known about the structure of CDM halos on the very small 100 pc) scales 
relevant to the formation of a spike[12]. 

2. Stellar Density Profiles 

The predictions of the adiabatic growth model can be tested against stellar luminosity 
profiles in galactic nuclei. The time scales for both physical collisions and gravitational 
near-encounters exceed a Hubble time for the stars in most nuclei[13]; in the absence 
of ongoing star formation, nuclei should have retained whatever density profiles were 
set up when the SBH gained its current size. Observed profiles are indeed well 
described as power laws within the SBH's sphere of influence, but with slopes that vary 
systematically as a function of galaxy luminosity[14, 15]. Faint ellipticals and bulges, 
My > —20, exhibit roughly the range of slopes predicted by the adiabatic growth model, 
1-5 < 7 < 2.5. However bright galaxies have shallower inner profiles, < 7 < 1.5. 

There are problems with identifying even the steep profiles in faint galaxies with 
the adiabatic growth picture[16]. Spikes produced by adiabatic growth almost always 
exhibit an inflection at r ~ (Figure 1); the logarithmic slope increases inwards. In 
real galaxies, the slope decreases inwards; inflections are virtually never seen, except 
at much smaller radii where they are associated with point-like or non-thermal nuclei. 
Avoiding an inflection requires fine-tuning of the initial conditions: either the initial 
profile has to be steep, 70 ~ 2, which obviates the need for adiabatic growth; or the 
SBH has to be just the right size that the pre-existing slope at r ~ r h matches the final 
slope 7. 

The rapid rotation, high densities and disky isophotes of low-luminosity ellipticals 
and bulges imply that gas dynamics played a dominant role in their formation [17]. The 
steep stellar density spike at the center of the Milky Way, p* ~ r~ 2 , may be the result 
of sustained star formation from a reservoir of dense molecular gas[18]. Semi-analytic 
models of galaxy formation[19] suggest that galaxies fainter than My ~ —20 were formed 
from gas-rich progenitors. These arguments suggest that the structure of nuclei in faint 
elliptical galaxies and bulges is a combined result of SBH growth and gas dynamics. For 
this reason, faint galaxies are not ideal testing grounds for the adiabatic growth model. 

A cleaner test of the model comes from bright elliptical galaxies, and here the 
model clearly fails: nuclear profiles in bright galaxies are almost always flatter than 
the minimum slope allowed by adiabatic growth, 7 ~ 1.5. Profiles are sometimes so 
flat that they are best described as "cores," regions of nearly-constant density; indeed 
central minima have been claimed in a handful of galaxies[20]. Core radii can be ~ 10 2 
pc or more. Something other than adiabatic SBH growth is required to explain this 
structure. 
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Figure 2. Mass deficit vs SBH mass for three different assumed values of 70, the 
assumed logarithmic slope of the density cusp before changes due to ejection of stars 
by a binary SBH. (a) 70 = 2; (b) 70 = 1.75; (c) 70 = 1.5. Solid lines are M def = M.; 
units are solar masses. (From Ref. 26.) 



3. The Binary Black Hole Model 

A natural way to explain the low central densities of bright galaxies is mergers. A SBH 
is effective at tidally disrupting the density cusp of an infalling galaxy [21, 22]. If the 
latter also contains a SBH, a binary SBH will form, which ejects passing stars or dark 
matter particles via the "gravitational slingshot." [23, 24] A binary SBH ejects of order 
its own mass before either coalescing in a burst of gravitational radiation, or, if the 
supply of ejectable matter runs out, stalling at a separation of ~ 0.01 — 1 pc[25]. 

The binary black hole model does a credible job of explaining the central density 
profiles of bright galaxies. Figure 2 plots the mass deficit[26]: 
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the difference in the integrated mass between the observed density profile, and a 
p*(r) oc r~ 70 profile extrapolated inward from the turnover radius 77,. There is a good 
correlation of Md e f with M., although the constant of proportionality can lie anywhere 
between ~ 1 and ~ 10 depending on the choice of 70. Larger ratios would be consistent 
with a picture in which galaxies form through a succession of mergers[25]. 

Interestingly, another class of stellar system exhibits both spikes and cores: the 
globular clusters. Globular clusters may contain "intermediate mass" (M. ~ 1O 3 M ) 
black holes, although the evidence [27] is not compelling. However they almost certainly 
contained an early population of ~ 1OM black holes, remnants of the first generation of 
stars [28], which would have spiralled to the center and displaced the lighter stars before 
ejecting themselves via the gravitational slingshot. This mechanism appears capable of 
creating the cores[29]. The steep central spikes seen in some globular clusters[30] are 
probably a consequence of collisional relaxation. 
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4. Implications for Dark Matter Spikes 

A dark matter spike could have survived at the center of the Milky Way if no mergers 
occurred since the era of SBH formation[31]. But the formation of the SBH was itself 
probably triggered by a major merger which channeled gas into the central regions [19]; 
indeed a binary SBH may have played a crucial role in the channeling [32]. It is hard to 
see how a dark matter spike could have survived such an event, and in fact the quasi- 
stationary conditions critical for the formation of the spike may never have existed[33]. 
These arguments suggest that a dark matter spike is unlikely at the center of the Galaxy, 
even if the Milky Way disk has not suffered a major merger in the last ~ 12 Gyr[34]. 
A conservative upper limit on the density of neutralinos at the Galactic center would 
be that given by CDM models of structure formation, ignoring the SBHs[22]. Spike- 
free density profiles might still allow testing of significant portions of MSSM parameter 
space [35]. 
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